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Abstract. We present a unified relativistic approach to inclusive electron scattering based on 
the relativistic Fermi gas model and on a phenomenological extension of it which accounts 
for the superscaling behaviour of (e, e') data. We present results in the A resonance region 
and in the highly inelastic domain and show some application to neutrino scattering. 



1 Introduction 

Electron scattering off complex nuclei is an ideal testing ground for modeling 
neutrino-nucleus cross sections, whose accurate prediction is necessary for the anal- 
ysis of on-going experimental studies of neutrino oscillations at GeV energies, 
presently being pursued in the MiniBooNE and K2K/T2K experiments [1,2]. 

Indeed inclusive electron (e, e') and charged-current (CC) neutrino {vi, l~) are 
closely related processes: in first Born approximation they involve the response of 
the nuclear system to a virtual boson, a photon or a W^, probing the electromag- 
netic and weak nuclear currents, respectively. 

Since the typical neutrino energies in oscillation experiments are of a few GeV, 
we will focus our attention on this kinematical domain. In this case the inclusive 
cross section shows a pronounced peak, the so-called quasielastic peak (QEP), 
at an energy transfer lo ^ \/ + mf^ — rriN, corresponding to the quasi-free in- 
teraction with the individual nucleons in the nucleus (here niN — nucleon mass). 
For high values of the momentum transfer q = \q\ and higher energy loss it is pos- 
sible to produce real pions and the cross section shows a second peak dominated 
by the resonant production of a Z\(1232) at w ~ \/ q'^ + rn\ — itt-n, where m/s. 
is the A mass. The width of these peaks is related to the Fermi momentum of the 
nucleons inside the nucleus and, in the case of the Z\-peak, also to the decay width 
of the A in nuclear matter Hence for a high enough value of q, these two peaks 
actually overlap and cannot be separated in inclusive experiments. At higher energy 
transfer the so-called second resonance region is found, where the iV*(1440)Pii 
(Roper), A^*(1520)£'i3 and Af*(1535)S'ii resonances are excited, evolving, at very 
high energies, into the Deep Inelastic Scattering (DIS) regime. 

Understanding the above spectrum in a unified framework has been the aim of 
some recent studies, which will be briefly summarized in this contribution. We shall 
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mainly address the electron scattering problem, where more data are available, and 
mention some applications to neutrino reactions. 



2 The A resonance region 

Since the kinematical domain we are exploring involves energy and momentum 

transfers of the order of (or higher than) the nucleon mass, the traditional non- 
relativistic approach is bound to fail and a relativistic approach to the problem is 
required. 

A fully relativistic treatment of the nuclear many-body system is a longstanding 
and extremely difficult task that is best pursued in special frameworks like the Rel- 
ativistic Fermi Gas (RFG) model, where basic symmetries Uke Lorentz covariance 
and gauge invariance are exactly respected. 

Hence we assume as a starting point an extension of the RFG model, which has 
been widely employed in the QEP, to the inelastic region. In this model, the virtual 
boson is absorbed by an on-shell nucleon described by a Dirac spinor it(h, Sh), with 
energy E^^ = i/h^Trn^. 

The basic ingredient of the calculation is the hadronic tensor W'"', which con- 
tains all the information on the nuclear structure and dynamics and yields the (e, e') 
differential cross section with respect to the lepton final energy e / and solid angle 
Qf according to the foUwing relation 

where a is the fine structure constant, e ; the lepton initial energy, rjfn, the leptonic 
tensor and = {lo, q) the four-momentum transfer. 

For the excitation of a stable resonance A''* of mass m* the RFG hadronic tensor 
turns out to be [3] 

where is the number of nucleons involved in the reaction, k = g/2mjv and rjp = 
V^pi^F + 2) = kp/mN stre the dimensionless transferred and Fermi momentum, 
respectively, and the "single-nucleon" tensor U*'^ embeds the information about the 
N -N* transition * . 

In analogy with the physics of the quasi-elastic peak [5], a scaling variable 
defined as follows 



ij)* = tp{K,, X;m*) = ± 




Actually in a fully relativistic framework the single nucleon physics cannot be exactly dis- 
entagled from the many-body part of the problem and the tensor Ul^" contains corrections 
due to the medium (see, e.g., ref. [4]). 
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with A = uj/2m]si and t — — X^, has been introduced. The quantity 

p{m*) = l + ^{m*^/ml-l) , (4) 

measures the inelasticity of the elementary process and reduces to 1 in the quasielas- 
tic limit m* = m^v- The physical meaning of the scaling variable -tp* is the follow- 
ing: ^fV'*^ represents the minimum kinetic energy required to transform a nucleon 
inside the nucleus into a resonance N* when hit by a photon of energy A and mo- 
mentum K. In terms of the scaling variable the RFG response region associated to a 
specific resonance is — 1 < V'* ^ 1 ™d the peak position corresponds to ij}* = 0. 

A realistic model for the resonance requires the inclusion of the decay width F 
in the hadronic tensor. This can be computed from the tensor W* ''(g, W) for a 
stable N* with mass by a convolution 

^ IIZ ^ (>v--..Tir(H^)V4 '^'''"-''^'"'»' • 

where the integration interval goes from threshold to the maximum value allowed 
in the Fermi gas model. 

The above expressions are independent of the specific transition under consid- 
eration: they can be used to describe the A, the iV*(1440) or any other resonance 
region [6,7] (after modifying to* and r{W) accordingly) as well as the quasielastic 
peak (in the F ^ Q limit setting to* ~ rnfj) [8]. 




Figure 1. Inclusive cross section per nucleon from ^^C. Llight blue: RFG without A width; 
green: RFG including the finite A width; blue: PWIA for the quasi-elastic peak; red: hybrid 
model obtained by adding the PWIA cross section in the quasi-elastic peak and RFG cross 
section in the A-peak. Experimental data are taken from [9]. 



As already mentioned the specific transition form factors are hidden in the tensor 
Ulf^ . Concerning the A resonance, a relativistic calculation of the (e, e') response 
in this region must use the full N ^ A vertex, which includes the magnetic (Ml), 
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electric (E2) and Coulomb (C2) excitation amplitudes. Although for years an im- 
portant program has been pursued to determine more accurately the quadrupole C2 
and E2 amplitudes in the Zi region (these being small compared with the dominant 
dipole Ml amplitude), our knowledge is still incomplete and has not been possible 
to undertake a full analysis in the sense of the work by Nozawa and Lee [10] of the 
effect of the C2 and E2 form factors in the nuclear zi-peak. Accordingly we use the 
parameterization [11] GEMiQ"^) = Ge.m{^)G'e{Q'^)I y'l - Q"^ /{i.hGeV/cf, 
namely we assume that the same dependence in is valid for the electric and 
magnetic form factors and that the isobar form factor falls off faster than the proton 
form factor. 

In Fig. 1 we show results for the nuclear inclusive cross section per nucleon from 
^^C compared with the experimental data. Light-blue hues correspond to the RFG 
for a stable A. The width is included in the green Une result according to eq. (5) 
and produces a broadening of the A peak and correspondingly a decrease of the 
strength. As an illustration of how one could improve the model in the quasi-elastic- 
peak region, we also show with blue Unes the quasi-elastic cross section computed 
with the semirelativistic PWIA model of ref. [12], which includes the momentum 
distribution of the finite-sized nucleus, producing the "tails" of the cross section, and 
the binding energy of the nucleons in the nucleus, yielding a shift to higher energies. 
Finally, we show with red lines the results computed with a hybrid model in which 
we add the PWIA cross section for the quasi-elastic contribution to the RFG result 
for the A contribution. 

As we can see in Fig. 1 , our results are below the data in the dip and A region. 
This was expected because other contributions coming mainly from two-nucleon 
emission and non-resonant pion production (not included in our model) also enter 
here. 

In Ref. [6] several relativistic effects and ingredients of the calculation were an- 
alyzed in detail. Here we just Uke to stress two basic findings. First, from a compar- 
ison between different Lagrangians in the treatment of the A excitation, it emerges 
that the Peccei Lagrangian [13], employed in the pioneering calculation by Mo- 
niz [14], is only appropriate for computing the transverse response for low momen- 
tum transfer, but in the longitudinal channel the fuU vertex [15] 

{A\j,\N) = ui (Cir^i^ + C2r',0 + c,rl^) u , (6) 

(where is the Rarita-Schwinger spinor describing a spin 3/2 particle) should be 
retained: indeed the Peccei Lagrangian, corresponding to the first term only, 
gives an unreasonably large longitudinal response. Second, we have found a large 
sensitivity of the longitudinal response to the inclusion of the Coulomb form factor 
of the isobar, especially for high q, a fact that could be of importance for investiga- 
tions of the longitudinal nuclear response. On the other hand, both longitudinal and 
transverse responses are found to be insensitive to the quadrupole E2 form factor. 

Let us now briefly comment on the second resonance region, namely where reso- 
nances heavier than the A are found. Among these, the A^* (1440)Pii (the so-called 
Roper resonance), occurring just above the Zi(1232), is particularly intersting, since 
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it can be viewed as a radial excitation of a three-quark nucleon state, analog of the 
breathing mode of the nucleus, hence carrying informations on the nucleon's com- 
pressibihty. A first step towards a theoretical description of the nuclear response 
functions in the Roper resonance region was made in ref . [7], where the correspond- 
ing RFG responses were calculated within several different theoretical models for 
the electroproduction amplitudes. In particular it was shown that, although the ex- 
perimental information was still insufficient to allow a stringent test of the various 
theories, the longitudinal response Rl associated with the Roper can be large com- 
pared with the contribution arising from the A near the hght-cone and that that the 
impact of the Roper on the Coulomb sum rule can be significant. 

3 The highly inelastic region 

To study what we call the "highly inelastic" domain, namely the region beyond the 
resonances, we assume that the final state can be described in terms of a recoiUng 
nuclear state plus a highly inelastic state of mass m*. In this case the RFG hadronic 
tensor can be written as [3] 

where pi and p2 are the kinematical boundaries at fixed q and w. Thus for each 
value of p (and hence m*) a "peak" can be identified, corresponding to the region 
— l<ip*<l, centered at t/j* = 0, whose width is a function that grows with 
q and decreases with to*. The single-nucleon inelastic hadronic tensor CZ/^g; can 
be parameterized in terms of two structure functions, wi and W2 (see ref. [3]). In 
computing it we employ phenomenological fits of the single-nucleon inelastic struc- 
ture functions measured in DIS experiments. Among the various parameterizations 
which can be found in the literature here we adopt the Bodek et al. fit of [17], which 
describes both the deep inelastic and resonance regions. 

Once the inelastic RFG modehng is in hand, we can go beyond this simple non- 
interacting model, still retaining its relativistic content. This we do by exploiting the 
so-called "superscaling" behavior of inclusive lepton-nucleus reactions. We shall 
see that this has a significant impact on the nuclear responses at high inelasticity. 

In order to introduce the concept of superscaling, we observe that Eq. (2) can be 
recast in the form 

WrflFG(9>'^>^*) = fRFGir) X G''-'{q,w,m*) , (8) 

where 

fRFGir) = lo{i-r'){i-r') (9) 

is the RFG superscaling function: hence by dividing the RFG nuclear tensor by 

an appropriate tensor G^'^ - whose definition follows immediately from Eqs. (2) 
and (8) - a universal function fupG is obtained, which does not depend on the 
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three variables q, u and kp independently, but only upon one specific combination, 
namely the scaling variable (3). 

As a consequence the RFG model predicts that if the inclusive cross section 
is divided by an appropriate function and plotted versus the corresponding scaling 
variable, no dependence on the momentum transfer q (scaling of the first kind) nor 
on the target nucleus, specified by the fermi momentum kp (scaling of the second 
kind), is found. The simultaneous occurrence of the two kinds of scaling is known 
as superscaling. 

In the quasielastic peak domain (m* = toat), superscaling has been widely 
tested against the (e, e') data in Ref. [18], where it was shown that both kinds of 
scaUng are fulfilled with a good degree of accuracy in the region < providing the 
momentum transfer is not too low. The scaUng analysis of the separated longitudinal 
and transverse responses has also proved that the scaling violations observed in the 
ip > region mainly reside in the transverse channel, whereas the longitudinal data 
do scale rather well in the whole QE region. 

However the shape of the experimental QE longitudinal superscaling function 
differs from the RFG parabolic one, extending outside the region — 1 < i/i < 1 and, 
most importantly, displaying an asymmetry around the QEP with a pronounced tail 
in the V' > region. The asymmetry of the superscaling function has been and still 
is the object of many investigations [19-23] and represents a stringent constraint 
that any realistic nuclear model aiming to describe (/, /') reactions should fulfill. 

An expression for a phenomenological QE longitudinal scaling function, 
JqeH'), was obtained by fitting the data [24]. Based on these results, we now make 
the following hypothesis: we assume that /q£;(V') provides a good description of 
/(^*) = = fritp*) ("scaling of the zeroth kind"), as it implicitly contains 

the initial-state physics, and thus we make, /or any m*, the following substitution: 

fRFoir) ^ fQEir) ■ (10) 

The corresponding results are illustrated in Fig. 2, where the inclusive cross section 
is shown together with the separated QE and inelastic contributions for a ^^C target 
at Einc = 4.045 GeV and 0e = 15° (a), 30° (b), 45° (c) and 74° (d). The calculation 
was performed in the Relativistic Fermi Gas including a phenomenological energy 
shift (red lines) and in the phenomenological extension of it (magenta lines), which 
we denote here as ERFG (Extended Relativistic Fermi Gas), based on the fit of the 
quasielastic scaling function Jqe [24]. 

We notice that for low scattering angle (a) the RFG model yields roughly the 
right position and height of the QE peak, but fails to reproduce the tails of the peak, 
giving in particular an unobserved dip at ~ 800 MeV. On the other hand the 
ERFG, while reproducing the data in the tails better, significantly underestimates 
the cross section at the peak. This is related to the fact that the peak of the phe- 
nomenological function /qe is lower than the corresponding RFG value. For higher 
angles [Figs. 2 (b),(c),(d)] the data lie roughly in between the predictions of ERFG 
(smaller) and RFG (larger) models, the former again reproducing the low-w behav- 
ior better. As a general result we observe that as the scattering angle increases the 
range of validity of the ERFG also increases. 
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Figure 2. Inclusive cross section for electron scattering from carbon at -Einc = 4.045 GeV 
and = 15° (a), 30° (b), 45° (c) and 74° (d), versus the energy transfer. Red: RFG; ma- 
genta: ERFG. The calculation includes an energy shift uisMft = 20 MeV and the separate 
QE and inelastic contributions to the cross section are shown (dashed). Data are from [16]. 



An important comment is in order. The RFG and ERFG models consider here 
the Ip-lh one-body contributions both for elastic scattering from a nucleon in the 
nucleus and for representations of the single-nucleon inelastic spectrum, thereby 
incorporating effects from meson production, excitation of baryon resonances (no- 
tably the A) and, at high excitation energies, DIS. However, in this region and be- 
yond effects arising from reaction mechanisms not included here, namely, those 
coming from correlations and both Ip-lh and 2p-2h meson-exchange currents can 
be also important [25-27] and, from some preliminary study, they tend to bring the 
total (the present ERFG contributions plus these additional MFC contributions) into 
better agreement with the data. Therefore the fact that the FRFG yields a cross sec- 
tion that is below the data is somehow encouraging, since this leaves room for the 
above-mentioned effects to provide the balance. 

Similar results are obtained for the superscaling function: in Fig. 3 the total scal- 
ing function / is shown as a function of the QE variable t/j' (the "prime" indicating 
the inclusion of a phenomenological energy shift, needed to reproduce the QEP po- 
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sition) for four different nuclei, within the RFG (left panel) and ERFG (right panel) 
models, at = 3.595 GeV and 6^ — 16°; experimental data are obtained from 
the measured inclusive cross sections divided by the single nucleon cross section 
and the curves are obtained by dividing the theoretical inclusive cross section by the 
same quantity. A closer inspection of the transverse superscaling function friip'), 
performed in ref. [3] at the same kinematics, shows that the discrepancy between 
"data" and "theory" is larger for the transverse case than for the total scaling func- 
tions at this scattering angle (6 = 16^). This indicates that extra contributions should 
be added to the nuclear model, going beyond the present one-body description, and 
that these must act mainly in the transverse channel. 
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Figure 3. Total superscaling functions /(V'') for Ee — 3.595 GeV and 9^ = 16°. Theoretical 
results obtained within the RFG are shown in panel (a), while the ERFG case is presented in 
panel (b). 



4 Scaling in the A region and application to neutrino scattering 

In the previous Section the superscaling function Jqe extracted from the quasielas- 
tic peak data has been used to calculate the full inelastic (e, e') spectrum. Two ba- 
sic assumptions underly this approach: first that the longitudinal and transverse su- 
perscaling functions coincide (0-th kind scaling) and second that the superscaling 
function, embodying the nuclear initial and final state interactions, is the same in all 
kinematic al domains, the latter being characterized only by the structure functions 
Wi and W2- 
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In Ref. [28] a different approach to the A resonance region was taken. First the 
contribution of the A has been isolated by subtracting from the total experimen- 
tal cross section the quasielastic contribution, reconstructed using the superscaling 
function /qe introduced above. Next the left-over cross section has been divided 
by the appropriate N A single-nucleon cross section and the result has been 
displayed versus the scaling variable tpA, given by Eq. (3) when m* = itia- The 
results are found to scale quite well [28] for tpA < 0, suggesting that this proce- 
dure has indeed identified the dominant contributions not only in the QE region, but 
also in the A region. Of course for tp/s. > higher resonances come into play and 
this procedure is no longer correct. The residual scaling function f'^, whose fit is 
plotted in Fig. 4, while similar to Jqe, differs in detail: it is somewhat lower, is 
shifted slightly and is more spread out over a wider range of scaling variable. This 
is not unexpected, since implicit in this approach is the fact that the A brings with 
it its own width and shift. A microscopic analysis of //i is presently being carried 
out [29] in order to deconvolute the width from the total response and see whether 
or not the underlying scaUng function is indeed the basic /qe deduced above. 




Figure 4. Fits of the superscaling function for the quasi-elastic (QE) and /^-resonance regions 
plotted versus the corresponding scahng variable t/j and compared to the RFG result. 



As mentioned in the introduction, a major advantage of the approach above il- 
lustrated is that the scaling functions /qe and //j extracted by (e, e') data can be 
used to predict neutrino-nucleus cross sections in both the quasielastic and A re- 
gions. Indeed, as tested in ref. [28] in a wide range of kinematical conditions, these 
phenomenological functions give by construction a good description of the elec- 
tron scattering data and are in this sense model-independent. The so-called SuSA 
(SUperScaling Approximation) approach amounts to multiply the two superscaling 
functions for the corresponding neutrino-nucleon elementary cross section. 

Several applications of the SuSA approximation to neutrino and antineutrino 
scattering are shown in Ref. [28]. In Fig. 5 we display, as an example, a comparison 
of the SuSA and RFG results, clearly showing that the relativistic Fermi gas badly 
overestimate the cross section in both the QE and A regions as compared to the 
phenomenological model based on superscaling. 
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Figure 5. Neutrino reaction cross section for = 1 GeV and 6 = 'X5 degrees. Red line: RFG result; 
black lines: SuSA result (dashed: A; dotdashed:quasielastic, solid:total) . 



It is worth mentioning that more fundamental approaches, based on relativistic 
mean field theory with final state interactions [20, 22] and on the coherent density 
fluctuation model (CDFM) [30,31] and aiming to justify the properties of the exper- 
imental superscaling function, have been recently carried out both in the quasielastic 
and in the A regions and applied to neutrino reactions, leading to results similar to 
the ones shown in Fig. 5. 

5 Conclusions 

We have shown how the full inelastic spectrum of inclusive electron-nucleus scat- 
tering can be described within in a fully relativistic unified formalism which can be 
applied in the few GeV energy domain. The model, based on the Relativistic Fermi 
Gas, takes into account initial and final state interactions in a phenomenological 
way through the use of a superscaling function directly extracted from (e, e') data. 
This allows one to make model independent predictions for neutrino-nucleus cross 
sections, needed in the analysis of v oscillation experiments. 
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